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Abstract 

We compute the renormalised energy momentum tensor of a free scalar field coupled to 
gravity on an n+l-dimensional Einstein Static Universe (ESU), Rx 5'", with arbitrary low energy 
effective operators (up to mass dimension n + 1). A generic class of regulators is used, together 
with the Abcl-Plana formula, leading to a manifestly regulator independent result. The general 
structure of the divergences is analysed to show that all the gravitational couplings (not just 
the cosmological constant) are renormalised for an arbitrary regulator. Various commonly used 
methods (damping function, point-splitting, momentum cut-off and zeta function) are shown 
to, effectively, belong to the given class. The final results depend strongly on the parity of n. A 
detailed analytical and numerical analysis is performed for the behaviours of the renormalised 
energy density and a quantity a which determines if the strong energy condition holds for the 
'quantum fluid'. We briefly discuss the quantum fluid back- reaction problem, via the higher 
dimensional Friedmann and Raychaudhuri equations, observe that equilibrium radii exist and 
unveil the possibility of a 'Casimir stabilisation of Einstein Static Universes'. 
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1 Introduction 

The Casimir effect [1] has been under intense experimental scrutiny over the last ten years, since the 
landmark experiment of Lamoreaux [2], and has been experimentally well established in different 
setups. A very interesting property of the effect is that it can be either attractive or repulsive, 
depending sensitively on the boundary conditions, geometry and topology of the background under 
study It would be desirable to have some physical understanding on what, exactly, are the 
conditions for the effect to be one or the other. In particular because a repulsive Casimir effect 
may have many interesting applications, as for example, to quantum levitation [Ij. However such 
understanding does not exist at present time. 

The issue of a repulsive Casimir effect is particularly timely in the context of gravitational 
physics, since repulsive gravity seems to be dominating our universe at present times [3 [6]. More- 
over, repulsive gravity is required to power an early inflationary epoch (see [7] for an example of 
Casimir induced inflation) and might play a role in avoiding the Big Bang singularity, through a 
bounce (see [8] for an example of a Casimir induced bounce). Thus, any quantum effects that might 
produce repulsive gravity are of interest, and the Casimir effect is an experimentally confirmed one. 
The effect has also great relevance in the context of Kaluza-Klein theories, where the compact- 
ness of the extra dimensions confines any quantum field, and hence produces a Casimir force that 
tends to make the extra dimensions perturbatively unstable (as was first argued by Appelquist and 
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Chodos [9]). In general it seems very hard to find Kaluza-Klein type configurations with compact 
dimensions which are stable under quantum effects (see [lOj . section 10), without invoking super- 
symmetry. By the same principle, if our universe has a compact topology, as in a A: = +1 FRW 
cosmological model or in A; = 0, —1 cosmological models with identifications, the Casimir effect will 
induce a force for any quantum field. 

Given the lack of general understanding concerning the sign of the Casimir force, one might 
consider families of models and examine the change in character of the Casimir force as one varies 
continuous or discrete parameters. Within this spirit we study the Casimir effect of a scalar field, 

with arbitrary mass /U and arbitrary coupling ^ to the Ricci scalar of the background on an 
n + 1-dimensional Einstein Static Universe (ESU) M x S*" of radius R. Thus, the model under stud 
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is a 4-parameter family, labelled by one discrete, n, and three continuous parameters, /x, ^, R 
We will present closed expressions for the renormalised energy density pren and pressure pren, as 
functions of the aforementioned parameters. These are obtained through a generic regularisation 
procedure, which can, effectively, be connected to several well known methods such as point split- 
ting, "momentum" cut-off, zeta function and different types of damping functions, including the 
one discussed in [H]. The mathematical manipulations are done by using the Abel-Plana formula. 

Let us note that special cases of our analysis have been performed in the literature. For instance, 
ref. [12] considers the n = 3 case and conformal coupling; ref. [13] considers general n, but only 
conformal coupling and zero mass. In both cases the results in these references agree with ours. In 
the first case, a self consistent Einstein Static Universe is found, taking into account the quantum 
effect; we find similar solutions (see section 4). In the second case, our computation reproduces the 
result of jl3] in the appropriate limit (see and discussion thereafter). 

Besides the intrinsic interest in understanding the behaviour of the Casimir force on a D- 
dimensional Einstein Static Universe we will be also interested in a particular physical application. 
It is well known that ESUs are, classically, perturbatively unstable towards radial perturbations. 
One may then ask if the extra Casimir force induced by the quantum field could give rise to stable 
equilibrium radii. To tackle this question one should use the semi-classical Einstein equations 

which yield the higher dimensional Friedmann and Raychaudhuri equations (see e.g. [TTj) 

i?2 + 1 2k^ 



R"^ n{n — 1) 
with the subscript 'cla' denoting the classical quantity, and 

R{t) 2 {ri - 



R{t) n{n - 1) 



"da ^ "ren 



(1) 



(2) 



where = 8ttG^^^^^ and 



ain) = pin) + ^^^(n) . (3) 

n — 2 



^ Other studies of the Casimir effect on n-dimensional spheres include (see also |10) . chapter 9) for a scalar 
field with n-dimensional spherical Dirichlet boundary conditions and [10] (chapter 10 and references therein) for a 
scalar field in a Kaluza-Klein setup of type iW* x 5*". 
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Thus, in this paper, we will focus our analysis on the renormalised quantities p^en and cTren, the latter 
determines if the Casimir effect gives an attractive or repulsive contribution to the gravitational 
force. In a follow up paper |19j we will examine in detail the dynamics of ([T|) and ([2]). 

This paper is organised as follows. In section 2 we present the theory under study, the unrenor- 
malised quantities and the renormalisation technique. We are then led to equations (|37p - (|40p . which 
are the main technical results in this paper. To gain insight about these very generic expressions 
we analyse their behaviour in section 3, by presenting numerical plots. In section 4 we present 
a summary of our results and observe that equilibrium radii exist, even including the quantum 
contribution; we also suggest the possible existence of stable equilibrium radii based on the results 
herein and ([1]) and ([2]), an issue to be studied more thoroughly in [19]. In the appendices we discuss 
some technicalities used in this paper and show how the general regularisation scheme used makes 
contact with commonly used regularisation methods (point-splitting, "momentum" cut-off and zeta 
function). 

2 Renormalised T^j, for a Scalar Field on ESUs 
2.1 General considerations 

We will deal with a simplified model containing the gravitational field, g^^, and a real scalar field. 
We assume that we are at sufficiently low energy so that these are the only dynamically relevant 
entities which emerge from a full theory, by integration of any other (more massive) degrees of 
freedom. Thus, the low energy effective action considered herein contains all possible combinations 
of these fields as allowed by the surviving symmetries. The coefficient of each term is unknown, 
a priory, and arises from the full theory. However we can still make some considerations about 
magnitudes by identifying the relevant scales in the full theory and in the problem at hand. 

The two scales to be considered are the Planck length (which is assumed to be the fundamental 
scale in the full theory) and the largest length scale (or lowest energy scale) in the limit we are 
studying. All other dimensionful parameters in the problem should be in between. We denote the 
former by Lpi and the latter by R. The most general Lorentz invariant effective action containing 
g^i, and ^ in (n+l) dimensions will be given by, 

S.,.. = I r^^xJ2a, (|)'^'^'^^!r/i...,0(;)-'"^ {d,^) , (4) 

where all lengths are measured in units of R (which is why they are denoted with a tilde) and 
the fields have also been redefined in units of the appropriate quantity constructed from R. C^^ -j 
denotes an arbitrary tensor operator constructed from $ and its derivatives, with (mass) dimension 
Qj, and TZ^P^^ is an arbitrary tensor constructed from the metric and its derivatives, with dimension 
Pj. The sum over j represents all possibilities for {pj,qj,kj}, i.e all possible Lorentz invariant 
combinations of T^-f^J^'s and O^'^ -j's. Finally Lj is a length scale for the order j operator, which 
must be restricted by Lp; < Lj < R and aj is a dimensionless coefficient. 

If all scales in the fundamental theory are of the same order as the Planck scale, then the L'^s are 
of order Lpi. However, more realistically we should expect other larger scales at smaller energies. 
These may be generated by the dynamics of the theory or some intrinsic symmetry. Thus, some of 
the Lj's may be of the same order as R. 
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To proceed we assume that only a finite number of terms in expansion ^ survive at low 
energies and that all higher order terms are suppressed. Therefore we must decide, in advance, a 
given order for truncation. Once we have chosen such an order, we must include all operators up 
to the corresponding mass dimension. All the coefficients are then low energy parameters to be 
fixed by experiment. Then we can go on to analyse the quantum dynamics of the fluctuations of 
the fields. When doing so, we have to keep in mind that j20] : 



1. The high energy contributions to any observable, must be suitably regularised in a way that 
does not affect the low energy physics. 

2. Any quantity which can be eliminated by a redefinition of the low energy parameters (ajL^^"*"* ) 
must be discarded. 

3. Any contribution which is of higher order in Lpi/R than the chosen order for truncation, 
must be discarded. 

2.2 The action 

Consider a scalar field with mass /i coupled only to gravity (i.e. with no self coupling or 
equivalently potential), through the Ricci scalar, TZ, of a curved space-time. The scalar field part 
of the theory is described by the action 



These are all the possible operators with mass dimension up to n + 1, which couple to gravity, 
contain no self coupling and that will contribute in our model. Indeed, the only remaining possible 
operators are of the form <I>(C -l-T^^^-')), where C is some constant. By a field redefinition these can 
be exchanged by operators of the form d^dTZ^^^^ and others which are purely gravitational. The 
latter may be absorbed in a redefinition of the gravitational couplings in ([6]); the former will not 
contribute to the scalar field equation of motion in our model, since all curvature invariants will be 
constants. We shall further assume that all other higher order operators are suppressed by powers 
of Lpi/R (which is equivalent to consider that all corresponding Lj's are of order Lpi). 

Analogously for pure gravity we must consider all gravitational operators up to that same order 
of truncation pj + qj = n + 1 (c.f. ([1])). Therefore the action for gravity is, 

. . / [(n+l)/2] \ 

where [. . . ] denotes integer part. Note that all gravitational operators are of even dimension. On 
dimensional grounds, from (jH), we expect all terms in ([6]) which are of quadratic (or higher) order 
in the curvature to be heavily suppressed by powers of the Planck mass (if we assume a single 
fundamental scale for gravity) H However, for consistency, we still have to include them in our 



^Moreover there are strong experimental constraints [18] on theories of gravity with higher order terms, which 
show that the extra coefficients must be small. We could also claim that there is some symmetry in the fundamental 
theory which sets them exactly to zero. Neither of these arguments is strictly necessary since we expect these terms 
to be suppressed in the first place. 
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effective theory, because they may get renormaUsed, even though we can neglect them for the 
purpose of looking at the evolution equations for the radius of the Universe at scales where the 
effective theory holds. 

Thus, the total action describing our effective theory is: 

'^Matter i (7) 

where the last term represents the contribution of the classical matter that will be considered to 
be a perfect fluid with positive energy density and zero pressure. In this simplified model we are 
neglecting the interaction of the scalar field with the classical matter. 



2.3 Bare energy density 

Now we look at the quantum fluctuations of the scalar field on the background metric. The equation 
of motion for <I>, obtained from ([5]), is 

(□ - $ = ^ 

and its energy momentum tensor 

T*, = a^$9,ci> + ^ (7^^^ _ D^D,) + g^, j^2e - [5,^'5"<I> + {fi^ + (TZ)^^] . (8) 

Conformal coupling is obtained by taking the coefficient = (n — l)/4n (and the theory is then 
conformal if /i = as can be checked by noting that T^^g^^^ = 0), whereas minimal coupling 
corresponds to ^ = 0. We will fix the space-time geometry to be the one of an ra + 1 dimensional 
Einstein Static Universe: 

ds%su^ = -df + R^dVLs^ , (9) 

whose Ricci scalar \sTZ = n{n — l) / B? , and where we have identified the radius of the universe with 
the scale R introduced just before equation This is a solution of d?]), if one takes $ = = a*^^^ 
and the cosmological constant and energy density of the perfect fluid to be related to R by 

2A n - 1 
n — z R'^ 

The eigenfunctions of the Klein-Gordon operator (□ — ^TZ — /U^) may be taken in the form 

ci>W =e-^-*y/_;;)(^) , (11) 

where the dots represent all necessary quantum numbers to completely characterise the hyper- 
spherical harmonic Y^^\6) on S"", with angular coordinates 6 (see, for instance [15tl23j). With this 
ansatz, we find the frequency spectrum 



^y£{e + n - 1) + ^n(n - 1) + 



' R 
with each frequency having degeneracy 



G No , (12) 



,(n) _ (2£ + n-l)r(n + ^-l) 

r(n)r(^ + i) 



= ^> . (13) 
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The unrenormalised vacuum energy density is obtained by a straightforward generalisation of the 
canonical quantisation procedure discussed in [8j; using the notation therein, we obtain the expected 
result: 



+00 



(«), ,{«) 



£=0 



where V^"^^ is the volume (strictly speaking, the hyper-area) of S"". 
For future notational convenience, we define at this point 



k 



R 



, T] = —5^ + ^n{n — 1) , 5 



n — 1 



thus p2p may be rewritten 



(n) 



Manipulating (|13p we can also rewrite it as 

n— 3 



(n) 



2kR 



(n-1) 



— TT 



i=0 



kR + 



n — 3 — 2i 



= P- (kR) , n > 3 



where we have introduced the useful notation 



P±ir) 



and the expansion 



(2p + 2)!J^^ 

- n 



(2p + 3) 



r p 



T2Q±(r2) , n = 2p + 3, 



TQ±(r2) , n = 2p + 4 



m,=0 
P+1 



m=0 



(14) 



(15) 
(16) 

(17) 



(18) 



(19) 



which defines the (positive) constants am and /3m- Note that the case ri = 2 is obtained with the 
product in (jl7p replaced by 1, whereas the case n = 1 is non-degenerate. 

2.4 A generic class of regularisations 

We now define a regularisation prescription for the energy density based on very weak and physically 
well motivated assumptions. This will allow us to separate and analyse the divergent and finite 
contributions. By using the first law of thermodynamics we can then obtain the regularised pressure 
and therefore the full regularised energy momentum tensor. In section \2M all the divergences and 
finite contributions which scale with the background radius R in the same way as the various 
gravitational terms will be subtracted, renormalising the corresponding couplings. 

To regularise the energy density we need to suppress the high energy modes. This is achieved 
by introducing a regulator, which should actually arise naturally from the fundamental theory by 
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integrating out the massive degrees of freedom; but our ignorance about such theory impUes that the 
way in which this regulator fundamentally arises is unknown. Besides suppressing the high energy 
(or momentum) contributions to physical observables, the regulator leads to a cutoff dependence in 
the bare couplings of our effective theory. But if this theory is to give a good low energy description, 
the infrared physics will be regulator independent, so such dependence is assumed to cancel when 
we renormalise the bare couplings 

We assume that, in the fundamental theory, there is a natural length scale 7-L which acts as 
regularisation parameter {L is a typical length scale in our model and 7 will turn out to be a very 
small dimensionless number). In addition we can introduce other dimensionless parameters, 7^, 
which could, for instance, control the form of the regulator before it is removed; alternatively they 
could be associated with ratios of other length scales in the problem. A realisation of the former 
possibility will be seen in appendix lA.2[ Thus, on dimensional grounds, a general regularised energy 
density which could arise from the fundamental theory is 

^ +00 

P^^^ = ^)Y.dt'-t\{lL.t\^.). (20) 
£=0 

The physical condition that the infrared physics should be regulator independent, amounts to say 
that the regularising function g, introduced in (pUj) . stays close to 1 for a wide range of energies 
around zero. This is equivalent to say that for sufficiently small energies 

— 2^6g <. — -fLujdog <. 1 , 
g to 

where 5 denotes a variation with energy and do a differentiation with respect to the first argument 
of g. Hence, to a very good approximation 

g{jLujf'\-fi) ~ g(0,7i) . 

Thus, an order 1 rescaling of 7 will not affect the low energy physics; this will turn out to be an 
important condition used in the isolation of the divergences, which is not surprising since 7 is to 
some extent arbitrary and any arbitrariness must be absorbed together with the divergences. 

In what follows, any remaining dimensionless parameters ji will not affect the reasoning so they 
are omitted; we therefore write 

+00 

P^^' = Ji^E4'^-^'''9{lL4'')- (21) 
£=0 

One way to interpret this formula is that only the energy spectrum is distorted. This is physi- 
cally reasonable since the background is taken to be fixed and thus the scalar field could still be 
expanded in hyper-spherical harmonics, so £ should still be a good quantum number. This type of 
regularisation is therefore consistent with spacetime symmetries. We will confirm it for the special 
case of point splitting regularisation in appendix lA.il 

It might seem that ()2ip simply introduces a damping function. It is certainly true that any 
damping function regularisation is of this form. But we will make the point that, if we leave the 
dimensionless function g unspecified - and we will obtain the renormalised energy density without 
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specifying g -, (|2T]) effectively includes a much broader class of regularisations. Indeed, as shown 
in appendix [X] for the special case of point-splitting, it will become clear that g does not even have 
to be a damping function; it might simply oscillate with increasing u;^ without ever decreasing its 
modulus. 

Since we expect 7L to be of the order of the Planck length and L is of the order of the typical 
length scales in our problem (which are either R or /i^^), 7 should indeed be very small, which 
justifies taking 7 to zero to identify the divergences. In fact, we will not have to set it to zero, since 
all cutoff dependent contributions will explicitly take a form which is readily absorbed in low energy 
couplings. The same will hold for all regularisations considered in appendix [XJ independently of 
any other extra parameters 7^ controlling the form of the regulator. 

For the calculations that follow, ()2ip is conveniently written in terms of a function /^(£ + 5): 

^ +00 +00 
2.5 Separating contributions with the Abel-Plana formula 

To separate the divergences from the finite part, we will use the Abel-Plana summation formul^ 

^/ ^ /'+°° ^, ^ , G'(O) . Giit) - G{-it) , 

Y^G{m)= G{t)dt + ^ + t Lt_\ dt, 

m=0 

in a way that generalises the procedure in |8]. The only restrictions on G{z) are that it must i) be 
analytical on the right complex semi-plane {z £ C : 5R(z) > 0) and ii) obey 

lim e-'^''y\G{x±iy)\ =0 . 

y-KX 

These are very weak requirements which are obeyed by most well behaved functions. Note that, 
when applied to (j22p , the second condition requires only that the function does not grow faster 
than exponentially as we go to infinity along the imaginary axis. This will enable us to obtain 
an answer which is independent of the regularisation and therefore consistent with any reasonable 
high energy theory which reduces to our effective theory. 

Applying the Abel-Plana formula to (j22p gives (with a simple manipulation that will become 
useful) 

MS) , , Mit + 5)-M-it + 5l^^ 



(23) 

where tmin = for > and tmin = \aR\ otherwise; I-^^div is a divergent contribution (in the 
7 — > limit) corresponding to the first integral and I-yjin is a finite contribution corresponding to 
the last three terms. Note that the renormalisation procedure will not simply amount to subtract 
I'y,div, c.f. section [2^ Using complex variable integration to simplify the latter (see appendix IB]) 



^For a thorough review of this formula and its applications see [21] 
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yields 



lim I-yJin 
7— >0 



dr 

aR exp(27rr) - 1 



ar- 



exp(27rr) + 1 



n = 2p + 3 



n = 2p + 4 



(24) 



for > and 



7-^0 ' ■' 



-IF 



+ 



exp(27rT) — 1 

^ Jo 



sin irt 
sin vrt 



cos vrt 



n = 2p + 3 



n = 2p + A 



(25) 



for < 0. To deal with the divergent contribution start by performing the change of integration 
variable 



u = RuJi 



in) 



to get 



-T duu'^^Ju^ -a?R'^Q-{u^ -a^R^)g(^v\ , n = 2p + 3 
2 J«o \ R J 



Uo 



duu^Q-{u^ - a^R^)g(^u] , n = 2p + 4 



where uq = aR for > and uq = otherwise. It is now convenient to separate the analysis 
according to the parity of n. 

• For n = 2p + 3 perform the change of variable 

u = \a R \v , 



to obtain 



4 p4 f+00 



L 



a^R 



■y,div 



dv - sQ-{s{aRf{v'^ - s))v'^g i-f\a\Lv) , 



where vq = 1 for > and vq = otherwise and s is the sign of a^. As argued in the last 
section, an order 1 rescaling of 7 should not affect the low energy physics; thus we have the 
freedom to perform the finite rescaling, 7 = |a|L7. Then we have 



4 7?4 r+00 



L 



a'^R 



dv \l — sQ-{s{aR)'^{v'^ — s))v^g (71;) . 



Thus, the divergent contribution in (|23p is obtained by putting all the pre-factors back and 
using the expansion ([T9|) : 

P+2 



E 

q=0 



A, 



R^i 



(26) 
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where 

A. 



m=m{q) 



2y, 



(n) 



i?=l 



p + 2-q 



/•+0O 

JVQ 



rj is defined by (fT5]l and m(g) = for g = p + l,p + 2 and m(g) = p — q otfierwise. It already 
becomes clear from (I26p that all divergent contributions scale with R in the same way as the 
gravitational operators in 

For n = 2p + 4 we simply perform the change of variable uL/R = v; extend the integration to 
V = and subtract the corresponding term, which exists only for > 0. Using the expansion 
for we get, 



P+i 

E 

q=0 



(27) 



where 



B. 



l)''/3p+i_g+fe fp+l-q + k 



k=0 



2V, 



(n) 



i?=l 



A: 



/• + 00 



and an extra finite term which is only present for > 0, 

Pextra — / ^ 



m=0 



(28) 



where we have defined the coefficients hm as 



^' (n) ^" 



(3 + 2m)!!T4'^\ 



and where is the volume of the n-sphere with unit radius. 



Assembling the results (l24|l . (j25[) - (l28[) . we finally get the full regularised result for (j23l) . up to 
order (7^^): 



rreg 



+ 



q=0 
p+1 



aR exp(27rT) - 1 



P+1 



hm{aR) 



2m+3 



k 9=0 



/22g 



m=0 



i?2p+5 



+ 



(_l)P+i f'R ^^P ^{T)^a'iR^ -t"^ 
exp(27rT) + 1 



RV^-^) 



dr 



n = 2p + 3 , 



, n = 2p + 4 , 



(29) 
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for a > and 



r p+2 

E 

q=0 



rreg 



|a|_R 



exp(27rT) — 1 
cos vrt 



- + 



+ / dtP_{ty-a?R'^-t'^ 



sin -nt 



5, 



E?Q 2RVi^) Jo 



'a\R 



dtP-{t)^-a?R^ -t"^ 



sin vrt 
cos Tit 



n = 2p + 3 , 



n = 2p + 4 , 



(30) 



for a? < 0. Now that we have the regularised energy density, the regularised pressure can be easily 
obtained via the first law of thermodynamics, as in [8j, 



p. 



(31) 



This is only a statement of conservation of energy (or equivalently covariant conservation of the 
energy momentum tensor) which is one of the good properties of our low energy theory that we 
want to keep. The result is 



pin) 
freg 



p+2 

E 

g=0 



n - 2(7)^,,^,;^,g,p , (-l)P P+{T){T^-a'R^ + ^l■'R^ 



n 



R'^1 



+ 



ni 



JaR Vr2 -a2i?2(exp(27rT) - 1) 



n = 2p + 3 



^ (n-2g)i?,,^,^,g,p ^ [{aRf^+ ^ - (2m + 3)(ai?)2-+ V'i?^] 



9=0 



m=0 



+ 



aR 



dT 



P+{T){a'R:'-fi'R'-T 



nRVi'^) Jo ^' V^?W^{exp{2TTT) + 1) 



n = 2p + 4 , 
(32) 



for > and 

P+2 

-E 

5=0 



„(-) 

Freg 



(n - 2g)^g,^,^,g,p ^ (-I)?' 



n 



R^i 



+ 



, P+(r)(r2-a2i?2+^2^2) 

«T — , ^ ; ^ - + 



Jo Vr2 -a2i?2(exp(27rT) - 1) 



+ 



1 



2nW(^ 



''l^, „ , ,/x2i22_a2^2_^2^^g^^ 



P+i 

E 

g=0 



(n - 2q)Bq^^^^^^. 



n 



R^i 



2nW(") Jo 



dtP^{t) 



V-a'^R'^ - ^2 sin vrt 

;u2i?2 _ ^2^2 _ ^2 gjj^^^ 

V-a2i22 _^cos7rt 



n = 2p + 4 , 
(33) 



for a2 < 0. 



2.6 Renormalising the energy-momentum tensor 

Now that we have separated all the divergences from the finite result, we proceed to renormalise 
the energy momentum tensor. To do it, any high energy contributions to the physical observables 
must be absorbed into the set of gravitational couplings in ([6|), which are fixed at low energies. 
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Such contributions are the divergences that we have identified in (p9|) - (|33|) : thus we need to show 
that they can be absorbed into the aforementioned gravitational couphngs. 
Varying the action ([Tj) we find the Einstein equations 



a' 

p=2 



where the only contribution to the scalar field energy momentum tensor is its vev {T^J). Due to 
the maximal symmetry of the spatial sections of the background Q, these equations reduce to a 
generalised Friedmann equation 

- A . + ' T "'"^ - [.S - ^fti] . (34) 

q=2 

where the explicit form of the functions F^'^\n, R) depends on the curvature invariants 

7^(29), and 

to a generalised Raychaudhuri equation 



The latter equation can be obtained from (j34p together with (jSip . Hence, to prove that we can 
renormalise all divergences it suffices to consider the Friedmann equation (j34p . It is then manifest 
that the scaling with R of the divergences shown in ()26p and (j27p is the correct one to be absorbed in 
a renormalisation of the cosmological constant {q = 0), Newton constant {q = 1) and gravitational 
couplings a^"^-* {q > 1). 

It is interesting to note that, by using a very generic regulator, we obtained almost all possible 
divergences. The single exception arises for even spheres for which the gravitational coupling of 
order q = 2p + 4 does not get renormalised (by a divergent contribution) . Which divergences show 
up depends on the specific regularisation method and hence on the high energy completion of the 
theory. 

Also interesting is the fact that we can accommodate theories with special symmetries which 
forbid some higher curvature terms in the gravitational Lagrangian, by choosing regulators which 
integrate some of the divergences to zero (or very close). Thus, this treatment is consistent, for 
example, with a theory of gravity which contains only the usual Einstein action. The fact that all 
the other couplings are small may be seen as an intrinsic property of the theory at the given energy 
scale. 

After this subtraction, we are left with a finite result. However, we must check if there are any 
remaining finite ambiguities. These would be additive contributions which are of the form i?^^*^ 
with a pre-factor independent of R. Since these could be absorbed in a renormalisation of the 
gravitational couplings they would lead to an ambiguity. The easiest way to look for such terms is 
to consider the asymptotic expansion of the energy density in the limit when i? — > oo. Equation 
(j3ip can then be used to obtain the asymptotic expansion for the pressure. With these expansions 
we easily identify the various subtractions that must be made, by looking at powers of R~^. 
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For the case of odd spheres (n = 2p + 3), the dominant contribution in the finite part of the 
energy density ([29j) is, in the large R hmit, 



(2p+3) 
reg.fin 



(/ii2)2P+5/2 



(36) 



27r(2p + 2)!iiy(") 

where '. . . ' denote sub-dominant terms; hence, for odd spheres, the dominating term is exponen- 
tially damped and no extra renormalisation is needed. The corresponding terms for even spheres 
(n = 2p + 4) come from the two finite terms in ()29p : 

P+l m+l 



(2p+4) ^ 
eg, fin 

m=0 fc=0 
p+1 



6™, 7/^^^™"^''+^ /^m + 3/2 



^2{p+fc-m+l) 



1) 



A;=0 



k 

1 

2 



22fc+2 



P. 



(2p+4) 
extra 



where were defined in (fT9]) : hence, for even spheres, we need to subtract p'ex-kra' order to get 
the final renormalised result. To obtain the term to be subtracted to the pressure we use ([3T|) . 
which yields 



„(2p+4) 
f extra 



p+1 m+l 

EE 

m=0 k=0 



m 



l)6™r?'=/i2™-2fc+3 + 3/2 
k 



{2p + 4)i?2(p+fc-m+l) 



A different way to argue that these subtractions are to be expected is as follows. On general 
grounds one expects the Casimir effect to be suppressed for very massive fields. Indeed, the mass 
introduces an effective infrared cut-off. Since the Casimir force results from the infrared behaviour 
of quantum fluctuations in some conflned setup, one expects the Casimir force to vanish as the mass 
increases. One can check from ()36p that this is indeed what happens for odd spheres. However 
the terms in p^extra^ yield a contribution for the energy that increases without bound for increasing 
mass. Similar terms are discussed in [3] (section 2.3 and 3.4) for the special case of the 2-sphere. 
Therein it is argued that an extra renormalisation must be imposed which eliminates them. In our 
setup such extra renormalisation appears naturally. 

Finally we can exhibit the renormalised energy density and pressure. For n odd they are 



^(2p+3) 
Pren 



-1)P 



dr 

aR exp(27rr) - 1 



(a^ > 0) 



UT- 



„(2p+3) 



exp(27rr) 



1 



+ 



1 



ni 



2i?y(") 

i?V^(") JaR ^' Vr2 -a2i?2(exp(27rT) - 1) 



MR 

/ dtP.{t)^-a'^R'^ -t"^ 
Jo 



cos TTt 

sin nt 



(a2 < 0) 



dr 



(a2 > 0) 



-IF 



— a'^R'^{exp{2TTT 



a'R' + ^2i?2^ 



+ 



1 



2nW(") 



-1) 
a^R^ 



-+ 



(37) 



(38) 



dtP-it y , 

^/^^^ 



COS vri 



^2 sin vrt 



(a2 < 0) 
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For even spheres the renormahsed energy density and pressure are 



^(2p+4) 
Pren 



p+1 

E 

m=0 



p+1 /.ai?, P^(^)^/^Tr2 



1 



a|_R 



dtP^{ty-a?R'^-t'^ 



exp(27rT) + 1 
sin vrt 



cos vrt 



„(2p+4) 
Pren 



g 6^ ((ai?)2™+3 - (2m + 3)(ai?)2'»+V2i?^ ^ 

m=0 



+ 



(_1)P+1 



n 







|a|_R 



ni?2p+5 
, P+(T)(a2i?2-^2^2_^2) 

or — = 

Va2i?2 _^2(exp(2vrr) + 1) 

;u2i?2 _ 0,2^2 _ ^2 gi^^^ 



dtP-{t) 



V-a2i?2 _ ^2 cosvrt ' 




„(2p+4) 
r extra ' 



(39) 



(a^ > 0) 
(a2 < 0) 



, „(2P+4) 
Pextra 



(40) 

Note that even though we have defined the renormahsed energy density using two different 
branches, the prescription for regularisation and subtraction of the divergences was exactly the 
same on both sides. Thus we have subtracted the same finite term for both branches in the renor- 
mahsed quantities. The fact that we have to divide in two branches comes only from mathematical 
manipulations; it is a consequence of our way of representing the finite contributions. 



3 Analysis of the renormalised vacuum energy momentum tensor 



Expressions (j37p -(|40p are, generically, non-trivial functions of one discrete and three continuous 
parameters {n, ^, p, R). In this section we will analyse in detail the behaviour of pi"n and aien 
-which play important roles in equations ([T]) and ([2|)- resorting to some analytical behaviours and 
numerical plots. 

A first observation to simplify the analysis is that the p dependence is quite simple. Indeed, it 
follows from ([371)- (001), together with that 

Pren and (Tren Can be written in the form 

(n) ^ ailUtpR) (n) ^ <^rel{^,pR) . . 

Pren ' ' 

indeed this follows from dimensional analysis. Changing variable to ^ = pR, we have 

(n) _ n+1 (^ren{^, R) (ji) _ n+1 '^ren 

Pren P Rn+1 ' ren P Rn+1 

Thus, we can simply consider the behaviour of pi'^n and Cren for ^ = which will be done 
separately for even and odd n in the next two subsections. The behaviour for any 7^ is simply 
obtained multiplying the vertical axis by a factor of /i""'"^ and dividing the horizontal axis by p. 

The massless case has to be considered separately, but the radial dependence is very simple in 
such case. One can use ([3]), (j3T]) and (0T]) to show that pi^n = pi^n/n and henc^ 

(n) ^ n - 1 

ren „ nrren ' \^^J 
n — 2 



^Note that this numerical value for fi has dimension of inverse length, since we are using units with h — c 
re-introducing these constants the numerical value of the mass is nh/c. 



^For n = 2 the role of a is played by a'^-* = 2p^'^\ which, in the massless case is cr'^' = p 



(2) 
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(n) (n) 

Thus cJren is simply proportional to pren- The latter function has a very simple R dependence; from 

P'^n = ^ . (43) 

Thus we just have to understand the behaviour of the function a.ren(0) which is displayed in figures 
[T]and[2]for the first few cases with n odd and n even respectively. Some generic features are: alen(0 
is always negative for small ^. For n odd aren(0 has (n — l)/2 zeros and goes to zero for large 



with sign (— For n even aren(0 has n/2 zeros and diverges - as can be checked from (|39 
- for large ^ with sign 

(since one of the zeros is an extreme) . 




Figure 1: alen versus ^ for /i = and n = 3 (one zero), ra = 5 (two zeros) and n = 7 (three zeros). 
Note that the sign oscillates with n for conformal coupling. 

A further special case is when ^ = = (ra — l)/4n (besides = 0), which corresponds to the 
conformal coupling. Then 

(-1)^' P+(r)r 

— r-r- / clT ; ^ , n = Zp + 6 

y Jo exp(27rT) - 1 ' ^ . (44) 



, n = 2p + 4 

The vacuum energy vanishes for even dimensional spheres. For odd dimensional spheres the vacuum 
energy approaches a finite non-zero value whose sign alternates with n, starting with a positive 
sign for (it is indeed negative for and zero for [3]). These behaviours for the conformal 
case were first observed in |13j . and can be confirmed in figures [U and [21 for the cases displayed 
therein. In particular, for the even n case, the zero at is the one corresponding to an extreme in 
each curve. 

Before proceeding to a numerical analysis of the behaviour of plen and ai^n for /i / 0, let us 
make a final observation concerning the p = case. Since the radial dependence of the renormalised 
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Figure 2: Oren versus ^ for = and n = 2 (one zero), n = 4 (two zeros) n = 6 (three zeros) and 
n = 8 (four zeros). Note that one of the zeros always corresponds to an extreme, which occurs for 
conformal couphng. 

energy density is now of the type l/i?""'"^ one could be led to think that such terms should be 
absorbed in the renormalisation of some gravitational coupling, following the discussion in section 
I2.6[ Note, however, that the renormalisation process was appropriately implemented considering 
^ 0, since we must consider all operators with mass dimension up to n + 1. In this case no further 
renormalisation was necessary. Once the final result was obtained one can study limits thereof (like 
the massless limit) without having to worry about further renormalisations. 

3.1 Odd n 

The first non-trivial case is n = 3, for which pren and cri^n are plotted in figure [3l as functions of 
R, and for a number of representative couplings ^. The behaviour of both functions is qualitatively 

(3) 

similar, and can be summarised as follows. For large radius preii goes to zero through positive 
values, in accordance with (j36p . For sufficiently small coupling one verifies the existence of one 
zero. The R value of the zero decreases with increasing coupling and, for sufficiently large coupling, 
the zero disappears and the curves become always positive. A consistency check is to verify ()44p : 
in this case, it implies that we must have prJn > for a = 0. Using (jlSp for this case we verify 
that a = corresponds to R = ^/l — 6.^. One can then verify that for all couplings represented, 

(3) 

Pren is indeed positive for the corresponding R. Note also that t he zero disappears before we get 
to conformal coupling. For coupling larger than conformal (.^ = 1/6 for n = 3), only the positive 
branches of (j37p and (j38p are needed to describe the behaviour of the renormalised quantities. 
Next we consider the case with n = 5, for which p^en and Uren are plotted in figures S] and 
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Figure 4: pi-Jn versus R for /i = 1 and C = 0; 5o> Tg' ^'^^^ ^"^^ similar figures, the rhs figure 

gives a detail of the Ihs one. 

El Again, qualitatively p^en and Uren have similar behaviours: for sufficiently small coupling there 
are now two zeros. Increasing the coupling, the two zeros disappear; first the one with positive 
first derivative and then the one with negative first derivative. The curves then become strictly 
negative. Again note that the asymptotic behaviour is in agreement with (j36p . To check (|44p . 

(5) 

which in this case implies that we must have pi-en < for a = 0, note that the two branches meet 
at ii = 2^/T^^^. It is simple to check that the renormalised energy density is indeed negative for 
all couplings represented at the corresponding radius. 

(7) (7) 

Moving on to the 7-sphere, one can see a pattern starting to emerge (for odd n). preii Sund (Jren 
are now plotted in figures [6] and [71 Again, the same (qualitative) behaviour is found for the two 
functions: one finds that, for sufficiently small coupling, there are now three zeros. As one increases 
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the coupling, the number of zeros decreases until none remains. The asymptotic behaviour is in 



which 



agreement with (j36|) . since pren approaches zero through positive values. To check 

(7) 

in this case implies that we must have pren > for a = 0, note that the two branches meet at 
R = Sy^l — 14^/3. It is simple to check that the renormalised energy density is indeed positive for 
all couplings represented at the corresponding radius. 

It seems tempting to generalise the observed pattern for any odd n sphere. One expects, for 
sufficiently small coupling, that pi^n and ai^l start off at negative values, possess (n — l)/2 zeros 
and vanish exponentially as — > +oo with sign 

(_l)(™-3)/2^ Continuously increasing the coupling, 
the two quantities start to lose their zeros, one at the time and alternate the sign for small values 
of R each time a zero is lost. For sufficiently large coupling (but always before conformal coupling) 
there are no zeros remaining and the sign becomes the same as the asymptotic one. 

(n) (n) 

A final observation relates the behaviour of pr-en and aren, as functions of R, for a given coupling 
in the massless and massive case: whenever the massive curve has an even (odd) number of zeros 
the sign of the massless curve coincides (is the opposite) with the asymptotic sign of the massive 
one. 



3.2 Even n 

We have already mentioned, at the beginning of this section, that for even spheres and a = = p 
the energy density vanishes. It is interesting to start analysing the leading behaviour in a for either 
branch of (j39p as a ^ 0. Denoting |a^i?^| = <C 1 in either branch we find (up to the term 
Pextra \ which is the same in both branches) 



„(2p+4) 
Pren 



{-If 
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(2p+4) 



R=l 



P0 3 

6 



32 



e^ + 0(e^ 



(a2 > 0) , 

(a2 < 0) . 
(45) 



For clarity we have displayed separately the expansions of the finite sum and the integral in (j39p . 
Remarkably, the finite sum - which appeared naturally in our renormalisation procedure - smooths 
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Figure 6: pi-Jn versus i? for /i = 1 and ^ = ^, ^, |, -g- 
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Figure 7: cJrIl versus for // = 1 and = 5j, ^, |, |- 



the energy density at a = 0. This is another piece of evidence for the correctness of our procedure. 

Let us start by considering the asymptotic behaviour of the energy density. An analysis of (j39p 
shows that the leading term as i? — > +oo is of order l/R'^^^^: 



(2p+4) R- 
rren 



-1)P 



4/xi?2y(2p+4) 



p+1 
m=0 



rj 



m+2 



2(m + l)(m + 2) 



+ (-irCfi(-2m-3) 1 



1 



22m+3 



(46) 

The second term inside the square brackets is always positive; the sign of the first term depends on 
rj. The simplest case is conformal coupling, for which rj = 0; then the sign of (j46p is (—1)^. This 
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will be readily checked in our numerical plots. For coupling smaller than the conformal value, rj 
becomes negative and the sign of (j46p might differ from (—1)^; indeed, we shall verify it is always 
negative for minimal coupling, regardless of p. 




Figure 8: piJn and alJn, = '^■Pren versus radius for /i = 1 and '^ = 0, 

The first non-trivial case is n = 2, for which p\.Jn and Uren are displayed in figure El A general 
observation is that, as for the odd n case, the qualitative behaviour of pren and cJren is always 
similar. For n = 2, the qualitative behaviour is also independent of the coupling: the curves are 
always negative, in agreement with (j46p . In particular there are no zeros. 




Figure 9: prJn and Oren versus radius for = 1 and '^ = 0,YgQ,^,^,^,^. 

The next case is n = 4, for which plen and firen are displayed in figure [9l For minimal coupling 
the curves are very similar to those for n = 2. But as one increases the coupling, the energy density 
and a become positive for sufficiently small radius, keeping their negative asymptotic value; thus. 
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within a range of couplings these functions have one zero, which is lost beyond a certain coupling, 
at which the asymptotic value changes sign. Indeed, specialising (j46p for p = we have 



.(4) R- 

rren 



96/xii2y(4) 



(8^ - 1) 



9 



+ 



449 
1680 



+ 



Thus, the asymptotic values changes sign for some ^ smaller than the conformal value, and beyond 
this coupling the curves are strictly positive. 




Figure 10: p^rJn versus radius for ^ = 1 and C = 0,^,^, ^, ^. 

The case n = 6 is considered in figures [TOl and [TTl For zero coupling, again, the curve resembles 
that obtained for n = 2,4. As one increases the coupling one observes the same behaviour discussed 
for n = 4, namely that the curve becomes positive, first for small R only and then for all R. But 
as one keeps increasing the coupling one now observes another inversion: the curves become again 
negative; first for small values of R only (hence possessing a zero) and then, for sufficiently large 
coupling, for all R. This is the asymptotic behaviour expected from (|46|) . 

The case n = 8, which is considered in figures [l2] and [13] confirms the emerging pattern. A 
further oscillation is then observed before conformal coupling, which makes the curves positive for 
sufficiently large coupling. Thus, for n = 2p + 4 one would expect the energy density and a to 
start off, at zero coupling, in qualitatively similar way to the cases analysed. As one increases the 
coupling, one will see p + 1 changes in sign. In each case, the sign change occurs first for small radii 
(hence allowing the existence of one zero) and propagates for all radii as the coupling is increased. 
After these p + 1 sign changes, the curves stabilise with sign (—1)^ for all radii, approaching zero 
as i? — > +CXD. Note that, in contrast with the odd n case, each curve can have, at most, one zero. 



4 Conclusions and Discussion 

Let us summarise the main results in this paper: 

1 - We have computed the renormalised energy momentum tensor of a scalar field with mass 
and coupling ^ to the scalar curvature on an n + 1 dimensional Einstein Static Universe with 
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Figure 11: Uren versus radius for // = 1 and ^ = 0) ^) ^) 




Figure 12: p^en versus radius for /i = 1 and ^ = 0)'§^i'^i^)'^)^;^)|§- 



radius i?. The resulting renormalised energy-momentum tensor is completely specified by the 
renormalised energy density and pressure, which are presented in (j37|) - (j4U| ). In the process 
we establish some technical results: 

1.1 - We show that different regularisation/renormalisation methods yield the same result, simply 

because the computation converges to a generic scheme, which effectively includes several 
regularisation methods commonly used, such as damping function, point-splitting, a "mo- 
mentum" cut-off and zeta function (c.f. section [2^ and the discussion in appendix E]) : 

1.2 - We show how, in these different renormalisation methods, various gravitational couplings 

(associated to gravitational operators with mass dimension up to n -|- 1) are renormalised 
and why all the renormalised couplings should be considered in a self-consistent effective field 
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Figure 13: ai^n versus radius for /z = 1 and ^ = 0)^)^)^)^!M'^'5i" 
theory treatment. 

We unveil a rather non-trivial, but clear, pattern for the behaviour of the vacuum energy density 
as a function of coupling, dimension, mass and radius. We also observe a transition, in the 
gravitational effect caused by the vacuum energy momentum, from repulsive to attractive. 
These features were studied in sectional wherein a detailed numerical analysis, also matched 
with some analytical behaviours, was performed for the renormalised quantities p^en and 0"^"™- 
These are the relevant quantities that enter the n + 1 dimensional cosmological equations ([T]) 
and ([2]). One can summarise the results of the analysis as follows: 

L - The generic radial dependence of pren and aien is given by (f4TI) . This renders the mass 



dependence trivial and one needs only to consider one representative of the massive case (say, 
^ = 1) and the massless case. In either case a pattern emerges from the analysis of the lower 
n cases and it is convenient to separate the analysis according to the parity of n. 

2 - Massive, odd n: The qualitative behaviour of pren and cjren is similar and displayed for 
n = 3, 5, 7 in figures [3]1TI For sufficiently small coupling, one observes that these quantities 
are negative for small R, possess (n — l)/2 zeros and vanish exponentially as R ^ +oo 
with sign (— 1)("~3)/2^ accordance with (j36p . Continuously increasing the coupling, the 
two quantities start to drop off their zeros, one at the time, and alternate their sign for 
small values of R each time a zero is lost. For sufficiently large coupling (but always before 
conformal coupling) no zeros remain and the sign becomes the same as the asymptotic one. 

3 - Massive, even n: The qualitative behaviour of pi"n and Cren is again similar, but distinct from 

the odd n case; these quantities are displayed for n = 2, 4, 6, 8 in figures [8lfT3l For sufficiently 
small coupling these curves are negative for all R, vanishing as a power law, l/R^'^'^, as 
R — > +00, in accordance with (j46p . Continuously increasing the coupling, n/2 — 1 changes 
in sign occur. In each case, the sign change occurs first for small radius (hence allowing the 



24 



existence of one zero) and propagates for all radii as the coupling is increased. After these 
n/2 — 1 sign changes, the curves stabilise with sign (—1)"/^ for all radii, approaching zero as 
R — > +00. Note that, in contrast with the odd n case, for each coupling there can be, at 
most, one zero for both p^en and dren- 

2.4 - Massless: The radial dependence becomes the simple power law ()43p . The only quantity to 

(n) 

be determined is the coefficient Oren, which depends on the coupling ^ and is displayed in 
figures [U for n = 3, 5,7 and figure ([2]) for n = 2,4,6,8. One observes that this coefficient is 
always negative for small ^. For n odd ai-ln{^) has (n — l)/2 zeros and goes to zero for large 
^ with sign (— For n even alen(0 ^/^ zeros and diverges - as can be checked 
from (|39p - for large ^ with sign (—1)'^/^. In this case, one of the zeros is an extreme. 

2.5 - Massless and conformal coupling ^ = ^con/ (conformal scalar field theory): the Casimir 
energy density vanishes for even n; for odd n its sign alternates as given by (I43p and (I44j) . 
This behaviour, although not the general analytic expressions, had already been noticed in 

m. 



We would like to close with a discussion of the backreaction of the vacuum energy density 
and stress and with the observation that a number of self-consistent ESUs can be found, in all 
dimensions, taking into account the quantum effects. 

Consider equations ([T]) and ^ with both classical and quantum terms, which are 

(") - A ^ = _L A. 

Pda - + j^n ' ^c/a n - 2 ' 

where A, ^ > represent the contribution of a positive cosmological constant and a pressure-less, 
positive density perfect fluid, and the quantum contribution is taken in the form ([^T]) . Equations 
([1]) and ([2]) imply that an equilibrium radius, Reg will be possible iff 

j^n-l ^ 2[Ao+^o + «^l(g,/io)] 

n{n - 1) (4^) 
2 A 
n 



where we have defined 



Aq — Ai?g ^-"^ , Aq = AReq , ^0 = fJ'Req 



To examine solutions of these equations start by considering the massless case, for which ()42p holds. 
Then, (|Tr|) becomes 

r: 



2 [Ao + (n + l)Ao] 



^"^ n(n - 1)2 

\ ' (48) 
= 2Ao - (n - 2)Ao 

n — 1 

One may think of these equations as determining, R(,q and ^ for given Aq^Kq^u. Generically there 
can be several solutions for each choice of Aq^Kq^u (or none). The solution can even be the same 
as the classical one pO|) . by fine tuning ^ to be one of the zeros of alen(C); such zeros exist for all 
n and, for n even, a special case is the conformal value. 
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For the generic case, (j47p . one may think of such equations as determining, R(,q and /io for given 
j4o, Ao,n,,^. Again, these equations can have a number of solutions, so that equihbrium radii can 
be found. An interesting question is if such equihbrium radii can be stable, in view of the classical 
instability of the ESU against radial perturbations. The quantum contribution seems, by itself, to 

(n) 

be able to have such effect, since the radial derivative at some of the zeros of cr^en is positive, thus 
producing a restoring force, by virtue of ([2]). But a dynamical analysis of ([T|) and ([2]) with the 
renormalised energy momentum computed herein must be performed to test this possibility. This 
is work in progress |19j . 

Acknowledgements 

We would like to thank Filipe Paccetti for very useful discussions. M.S. is supported by the FCT 
grant SFRH/BD/23052/2005. R.R is supported by the FCT grant SFRH/BD/35984/2007. This 
work was also supported by the FCT grant POCTI/FNU/50161/2003. Centro de Fisica do Porto 
is partially funded by FCT through POCTI programme. 

A Contact with commonly found regular isat ion procedures 

We have already seen in section [27il that the general regularisation scheme used, expressed by (|2ip . 
makes close contact with a regularisation scheme that employs a damping function. We now make 
contact with other well known methods commonly found in the literature. We will show how, using 
these methods, at a certain point of the computation one arrives at a special case of the general 
regularisation procedure considered in (j2ip . The main point is that all of them seem to collapse 
into the same result, once we have introduced a sensible cut-off scale. 

A.l Point-splitting, an alternative formulation 

Point-splitting constitutes a covariant regularisation scheme, thus naturally respecting spacetime 
symmetries [181 122] . The conventional approach to this method consists in regularising the energy 
momentum tensor by noticing that the divergences arise from evaluating operators (which are 
quadratic in the field) at the same point. Therefore, one simply evaluates such operators at 
different points. Using ([8]) we can formally define, 



where x is replaced hy x' ^ x in all primed quantities. 

To make contact with the treatment in section [2741 the first step is to identify the regularisation 
length scale 7-L. This issue is intrinsically related to the issue that equation (j49p is only valid as long 
as we refer to a particular coordinate system and work in components; it has no frame independent 
meaning - it is only a tensor in the limit of coincidence^ However, as argued in section 12.41 we 

®In the conventional approach, one writes the point splitted energy momentum tensor in terms of Green's functions 



hm (1 - 20 d^^d,.^' + i ($'7^^.$ - ^'D^D,^ 




(49) 
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expect the fundamental theory to have an intrinsic regularisation which produces a regularised 
energy momentum tensor. 

To construct a true point splitted regularised tensor we proceed as follows. We start by consid- 
ering x' as a fixed map, which plays the role of a regularisation parameter. To avoid confusion we 
define such a map as 

h: M -^M 

X h{x) , 

where M is the background manifold and a; is a point on the manifold. This map induces a pullback 
map, which associates a tensor at x to every tensor at h{x). By taking the map h to be the 
identity we recover the original tensor field. Therefore this is a very natural way of constructing a 
point splitted field. 

For a scalar field, this procedure is particularly simple since the pullback is given by composition 
of the scalar field with the map 

^h- M ^ M 

X ^h{x) = $(/i(x)) . 

So if we go back to ([8]) we can obtain a true point splitted energy momentum tensor by replacing one 
of the <I>(a;)'s by ^h{x) in each quadratic term in <I>. Similarly we could have started by performing 
the substitution at the level of the Lagrangian and consider it the regularised Lagrangian for the 
theory. Functional differentiation would then yield the energy momentum tensor. In any case the 
result is exactly the same 

T*/ = (1 - 20d^^d,^h + i {^hn^u^ - ^hD^D,^ - ^D^D.^h) + 



+ g^,{2i-^[d^^hd''^ + ^h{n^ + in)^] , (50) 



with the bare energy momentum tensor being recovered by taking 

T* = lim T*''' 

where id is the identity map. 

Expression (|50p is a manifestly covariant quantity and can be put in a form which makes contact 
with the conventional point splitting technique, if we express the covariant derivative terms in new 
coordinates x' = h{x), i.e. 

we obtain an expression similar to (|49p 



(51) 



(and its derivatives) for the scalar field. The objects arising can be given a meaning as bitensors and may be 
manipulated in a way to separate various contributions which transform as tensors at x (they also contain the 
geodesic distance between x and x' , for more details see ^22; ) . 
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where = ^{x') = ^h(x). The extra Jacobian factors make this quantity a true tensor; so in fact, 
we have shown that the point-spHtting technique can be put in a truly frame independent way. 
The conventional point splitted energy momentum tensor (j49p is obtained if 

thus, the conventional procedure just amounts to take the map h to be an (arbitrary) constant 
translation. The regularisation length scale 7L can now be seen as the proper size of this constant 
translation vector. 

We can proceed to compute the vev of T^^, by applying the standard procedure, as described 
in [18] , and make contact with the generic regularisation in section 12.41 The usual method consists 
of expressing such a vev in terms of the Hadamard elementary function. Here, however, we prefer 
to use an expansion in mode functions together with the action of the differential operators on the 
latter. For concreteness we focus our attention on the 3-sphere, for we can simplify the treatment 
using group theory. Write the metric as 

where {^r) c'-re left (right) forms on SU{2) ~ 5^0 The isometry group is six dimensional; the 
Killing vectors form the Lie algebra of 5*0(4), which decomposes into su(2)i x su(2)^. Thus we 
can always choose mode functions transforming under irreducible representations of the algebra. 

For the 4-dimensional ESU, consider the (complexified) Killing vectors {ko, k^, k^} where kg = 
—idt is the timelike Killing vector and the remaining k's are the dual vector fields to the ex's. The 
full isometry group of the ESU has Lie algebra £ = u(l) x su(2)i x su(2)^. Thus, functions which 
transform as irreducible representations of £, can be chosen as mode functions. Note also that the 
Casimir operator k^ of each of the su(2) algebras is in fact the same, so the action on a generic 
basis function \LO,£,mL,mji) = \uj) \£,mL,mR) is 



ko\uJ,i,mL,mii) = u \uj,e,mL,mR) , 

\uJ,i,mL,mji) = C+{i,mL)\i^,i,mL + l,rnR) + C-.{i,mL)\oJ,i,rnL - l,rnR) , 
k2 \uj,i,mL,mR) = -iC+{£,mL)\ojJ,rnL + l,rnR)+iC-.{£,rnL)\oj,i,mL-l,'mR), 
k^ \uj,£,mL,mR) = mL\uj,£,mL,mR) , 

k^ \uj,£,mL,mR) = £{£ + l)\uj,£,mL,mR) , (52) 
with similar expression for the R algebra (interchanging L and R) and where 



C±{£,m) = ^y£{£+l) -m{m±l) . 

The quantities {lo £ M,£ £ Nq} classify the sub-algebras and are interpreted respectively as 
energy and orbital angular momentum eigenvalues and mR^L G {~^)~^ + Ij--?^} are the re- 
spective azimuthal quantum numbers. Note that = e~*^* and that the spatial functions 
^aix) = \£,nT-LT'iTT'R) ohej the same differential equation as the hyper-spherical harmonics. The 



'^To see an explicit representation in terms of Euler angles see [25] , appendix A. In the notation therein, the 
complexified Killing vectors used herein are ]if = —i^f and = i^f. 
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index a denotes the set {i,mL,m]i\ and x is a unit vector in the embedding space of S"^ (1^^), 
defined by the angles that we choose to parameterise the sphere. 

By taking into account the dispersion relation, we finally recover the original eigenfunctions of 
the Klein-Gordon equation 

1 



(53) 



where we have assumed normalisation of Ya{x) on the 3-sphere, 

j d^xY^{x)Y*,{x) = 6aa' , 
and of <pa{t,x) with respect to the Klein-Gordon scalar product, 

(fix {(t)a{t,x)ko(f)*c,,{t,x) - 4>a,{t,x)ko(l)a{t,x)} = 6^ 



The advantage of choosing a non-coordinate basis of Killing vectors is that all the partial derivatives 
are replaced by Killing vectors which have a simple action on the basis functions. The point splitted 
energy density then takes the form 



2r' 



<I>,h 
ab 



-(1 - 20ka#kfe/$' + ^ (7^a6$$' + ^K>kb>^') + 



9ab (2^-2 



/"(5f kc$kd'«>' + (^^ + iU)^^' +x^x' , (54) 



where ka denotes the set {ko,k^} or {ko,k^} and x ^ x' denotes symmetrisation with respect 
to interchange of primed with unprimed variables (including indices). Note that we have denoted 
some of the a and h indices in primed form to avoid the use of Kronecker deltas which would make 
the notation denser. Finally, after some algebra, the vev of the energy density in the vacuum follows 
from using ([52]) 

Ph= 2'^^l{<t^a{t^x)4)l,{t' ,£)+X ^ x] . 

a 

This result is very simple and suggestive. It is simply a sum over all modes of the frequency squared 
times the (symmetrised) point splitted amplitude for the respective mode. Note that, since this 
statement does not depend on the number of dimensions or the basis of functions, it should be 
easily generalisable for all spheres. Note that ()53p also holds for any sphere, so in addition the sum 
over i can be factored out to get 



where 



(n) 
Ph 



ge{x,x') = cos [uj^^\t - t')) 2-F1 



^ +00 

yr^J2df^uj'l'^ge{x,x') , 



2y(") 



(55) 



. . ,11 X ,X 

-i,i + n-l,i+-- ^ 



we have used the generalised addition theorem for spherical harmonics 



E yc.{x)Y:ix') 

a\{e} 



2£ + n-l 
(n-l)yW^ 



(x.x 
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and a relation between the Gegenbauer functions (C^) and the hyper-geometric function 2F1 
Ctiz) = f-., . + 2A, A + V 



T{2\)T{v + lY V ' 2' 2 

where the only restriction is that A is not a negative semi- integer. Equation (j55p says that the 
regularised energy density falls into a similar type to the generic regularisation considered in the 
text (j2ip (note that when x — > x' ). The Abel-Plana formula can too in this case be applied. 

The quantity obtained remains finite when x = x' which shows that a timelike point separation 
still regularises the result. Since any finite contribution to the final result cannot depend on the 
map we are free to make such a choice. Thus, finally evaluating 2-^1 at x = x', which is 1, we 
obtain that (j55|) is exactly of the type of ([2T|) . with 



gi'yLiv'f''') = cos (^7^"^!"^ 

and therefore we have shown how the point splitting technique overlaps with our generic class of 
regularisations. All the results in section [2] follow straightforwardly. 

A. 2 "Momentum" cut-ofF and a class of damping functions 

A more explicit regularisation scheme would be to apply a frequency cutoff, which is in the same 
spirit as the standard momentum cut-off used in QFT. It would consist on stopping the sum at a 
given i = N which is then taken to infinity. This still contains the physical condition that high 
frequencies are suppressed, even though such a rigid regularisation seems unlikely to occur from 

(n) 

the full theory. The cutoff in induces a cutoff energy ojj^ which we identify with the inverse 



scale 7L and can be implemented by introducing a step function, 9{z), in (pTjl 

The application of the Abel-Plana formula require using an analytic continuation of the function g; 
thus we should take the Heaviside step function as a limit of a family of smooth functions. Taking 

5(7^^^^) = (n),,^ > (56) 

1 _|_ e~(^~'^« V^AT )/* 

for fixed t, we have a family of damping functions. The rigid momentum cut-off is obtained in the 
limit t — > 0. The results in section 2 now follow for all t. 

Alternatively, we could use the following generalisation of the Abel-Plana formula (see e.g. [21] ) 

N ,N G(Ar) ^ 



V G{m) = / G{t)dt 



^ G(0) ^ , Gjit) - G{-it) + G{N + it) - G{N - it) ^ ^^^^ 
2 In e^'^* — 1 







In principle we could do similar manipulations as before and obtain analogous identifications for 
the divergences when applying ([57|) to (fT4|) with the sum up to N. 
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A. 3 Zeta function 



Zeta function regularisation of (|14p is similar in spirit to dimensional regularisation, since it relies 
on analytic continuation. It is implemented by writing 



where C represents the appropriate zeta function describing the spectrum and q some constant 
factors. Commonly found cases are the Riemann zeta function, the Hurwitz zeta function and the 
Epstein-Hurwitz zeta function. The renormalised energy density is then obtained taking the limit 



and using the analytic continuation of the appropriate zeta function. In many cases the result is 
immediately finite; the analytic continuation of the zeta function "swallows" the divergences. But 
it might happen that the analytic continuation of the zeta function has poles at s = —1, in which 
case one needs to discuss the surviving divergences. As an example we consider below n odd and 
generalise the method applied in [8], which uses the Epstein-Hurwitz zeta function, to treat the 
case with > 0. 

The regularised expression (jSSp can, alternatively, be treated as a special case of the general 
treatment in section 2.4. To do that with the single parameter 7 used therein we relate the 
parameters u and s by taking s + 1 = ijf"^ jv. Then, we note that (|58p can be written in the form 



where 7L is identified with v . Thus removing the regulator s ^ — 1 corresponds to ^ +00 and 
thus to 7 ^ 0, as expected. The function g becomes a damping function to which the treatment 
in section 2 can be readily applied; the results therein follow. 

A. 3.1 Standard analytic continuation 

The typical procedure in zeta function or dimensional) regularisation is analytical continuation. 
The latter often hides the divergences and their physical interpretation as we have emphasised. For 
completeness we illustrate how, for n odd and > we can perform the manipulations which lead 
to the analytical continuation and check that we recover the same result, albeit represented in a 
different way. 

Using ([16]) and p7|) . ([58|) becomes, after a simple manipulation, 




(58) 



Pren 



s 




(f2Tl) . with 





which can be written as a linear combination of Epstein-Hurwitz zeta functions 




(59) 
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plus a finite sum 

„(2p+3) 



s+1 



■ p+1 



.m=0 



n=l 



X] CmCEH[i: + m-p-l, {aRf] + 



(60) 



P+i 



.m=0 



wliere tlie coefficients Cm are given by 



Co 



1 , 



(-ir 



+ (aRY 



m G N . 



ri,r2, -- ,rm.—0 

ri<r2<---<rm 



One can show that the last term in (j60p vanishes, in the limit s — > — 1. As long as > we may 
use the analytic continuation of the Epstein-Hurwitz zeta function (see e.g. [JJJ), obtaining the 
result 
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Hren 
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' 2m— 71— 1 ^ 



r(^^) 



(ai?) 



-2m+n+l 



for n ^ 3 and > 0, where Kq, are modified Bessel functions of a*'* order. The last term in ([6T 



is divergent, due to the poles of the gamma function. It should be dropped out and understood 
as a renormalisation of the gravitational couplings of all the terms in ([6j), except the cosmological 
constant. This illustrates how different renormalisation methods give rise to different renormalisa- 
tion of physical quantities, since with a damping function method of the type used in |14t [S] only 
the cosmological constant gets renormalised (a result that generalises for arbitrary odd n). This is 
a consequence of the different UV completions of the theory obtained from different regularisation 
methods. But for a renormalisable theory the finite contribution should be independent of the 
regularisation method, which is what we are obtaining here. 

Subtracting these infinite contributions one obtains the renormalised energy density. In the 
special case of the 3-sphere it reduces to the result obtained in [8], where it was shown that it 
coincides with the result obtained with a damping function method. That (j6ip yields, in general, 
the same result as the one obtained before with our more generic regularisation scheme can now 
be checked numerically, case by case, which we have done for the 5-sphere and 7-sphere. 

Finally let us mention that, for ^ = and conformal coupling, we have checked that the analytic 
continuation of the Hurwitz zeta function yields 



B Simplifying the integral I 

In this appendix we simplify the integral 

f^{it + 5)- f^{-it + 6) 



I 



o2irt _ 2 



dt , 
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Figure 14: Integration contours in (j62p . C1/C2 is along the positive/negative imaginary axis. 

which appears in (1230 . leading to the expressions in p4p and (I25p . First consider the complex 
variable z = it. The integral / becomes the contour integral 



where the paths of integration are exhibited in figure [TH We will deform these contours of inte- 
gration into closed contours enclosing regions wherein the integrands are analytic functions. An 
application of Cauchy's theorem then yields the integral / in terms of the integrals along the re- 
maining paths of the closed contours providing a manifestly real expression for /. To do this rewrite 
fj given in (p2|) as 



where P_ [t) is given by (jlSp . It is convenient to analyse separately the cases with non-negative 
and negative respectively. Note that we will implicitly use the fact that the regulator is analytical. 
For convenience, we will keep 7 7^ for some integrals; this is allowed for all the extra contributions 
are of order 7*^ with k positive. 

-k > 0: The denominators in the integrands of (I62p have zeros at z G Z. These are all compen- 
sated by the zeros of the numerators, except for the one at z = 0. Therefore the integrands 
have a simple pole at z = 0. Moreover, we have branch cuts as z = —5 zizir, t £ [\clR\, +oo[, 
which are induced by the branch cut of ^/z at z G M^. Thus, in this case we consider the 
contours of integration exhibited in figure [TS] (left). Let us consider separately the several 
contributions: 

Paths Ai and A2: Integrate along z = ziziA + x respectively, with x s] — 5 + e, 0]. The 
exponential factor in the denominator will suppress these contributions when A — > +00. 




(62) 




33 





Figure 15: Closed contours of integration: left /right figures correspond to the cases with > 
and < 0, respectively. 



Thus 



lim {Ia^+Ia^) = 0. 

A— »+oo 



Paths Bi and B2: Take A +00 and integrate along z = ±iT — 6 + e, with r G [aR, +00 [, 
respectively and e infinitesimal; taking the limit e — > one finds that the sum of these 
contributions is non- vanishing only for n odd, 



hm + Ib2 



!/ 

< JaR 

. 0, 



~ P-(zr)Vr^-Q^i?^ 

ar — ; : , n = 2p + 6 , 



exp(27rT) — 1 



n = 2p + 4 



Paths Di and D2: Integrate along z = zizir — 6, with r € [0, aR] respectively. One finds that 
the sum of these contributions in non-vanishing only for n even, 



Idi + = ' 



0, 



n = 2p + 3 , 



['I 

K Jo 



^-^^^ — r, / , n = 2p + 4. 

exp(27rr) + 1 



Paths El and E2: Integrate along z = —t, with t g]0, 5[, 

I El + 



f dyf^{y) 
Jo 
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Pole contribution at the origin: Integrate along z = ee^*^ respectively, with G [7r,7r/2] and 
e infinitesimal, to obtain for the sum of these contributions is, 

lim (/p,+/pJ = M) . 

Applying now Cauchy's theorem and substituting in ()23p one obtains the manifestly real 
expression for the renormalised energy density ()24|) . 



★ < 0: The pole is the same as in the case > 0, but the branch cuts are different. They are 
exhibited in figure [15] (right). Let us consider separately the several contributions: 

Paths Ai and A2: Same as in the case > 0. 

Paths Bi and B2: Take A +00 and integrate along z = itir — 6 + e, with r S [0, +oo[, 
respectively and e infinitesimal; taking the limit e — > one finds that the sum of these 
contributions is non-vanishing only for n odd. 



lim + Ib2 



e-*0 



dr — ; : , n = 2p + 3 , 



exp(27rr) — 1 



, 



n = 2p + A. 



Paths Di and D2: Integrate along z = ibie — 5 + t, with t G [0, \ a\R] respectively. In the limit 
e ^ one finds that the sum of these contributions is 



lim + Id2 



( 1 /"l*^!^ , COSTTf 

-- / dtP_(t)x/-a2i?2_t2!:^ ^ n = 2p + 3, 
2 Jo smvrt 

1 MR 



K 2 



dtP^{t)^J-a^m -f^ 



sin Tit 
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n = 2p + 4 . 



Paths El and E2: Integrate along z = —5 + 1, with t £ [|o|-R, 6[. One finds a similar result as 
in the case > 0, 



I El + 



aR\ 



dyfjiy) 



Pole contribution at the origin: This is the same as for a'^ > 0. 

Applying now Cauchy's theorem and substituting in (j23p one obtains the manifestly real 
expression 
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